First search for gravitational wave bursts with a network of detectors 
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We report the initial results from a search for bursts of gravitational radiation by a network of five 
cryogenic resonant detectors during 1997 and 1998. This is the first significant search with more than 
two detectors observing simultaneously. No gravitational wave burst was detected. The false alarm 
rate was lower than 1 per 10 4 years when three or more detectors were operating simultaneously. 
The typical threshold was H ~ 4 X 10 -21 Hz" 1 on the Fourier component at ~ 10 3 Hz of the 
gravitational wave strain amplitude. New upper limits for amplitude and rate of g.w. bursts have 
been set. 



The direct detection of gravitational waves will be a 
watershed event for both the physics of gravitation and 
the investigation of compact astronomical objects. A 
variety of astrophysical events are expected to produce 
gravitational waves of short duration (<C 1 sec), or gw 
bursts, such as the gravitational collapse of stars or the 
final few orbits and the subsequent coalescence of a close 
binary system of neutron stars or black holes . Due to 
the inherent weakness of such signals, and the difficulty 
in distinguishing them from a myriad of noise sources, 
the direct detection of a gw burst will require coincident 
detection by multiple detectors with uncorrelated noise. 
Searches for gw bursts over periods of observation of 1 — 3 
months have been performed in the past by pairs of cryo- 
genic resonant bar detectors [§-0] , setting upper limits on 
the incoming rate. A few days of observation have been 
reported for simultaneous operation of three cryogenic 
bar detectors || and, with much less sensitivity, of a pair 
of short-arm interferometric detectors ||. Upper limits 
on gw signals from coalescing binaries has been recently 
reported also by a single interferometric detector for 25 



hours of observation || . 

In the last few years, the increase of the number 
of cryogenic resonant detectors in simultaneous opera- 
tion has greatly improved the prospects of obtaining a 
confident detection of gw bursts. There are now five 
operational cryogenic bar detectors: ALLEGRO (Ba- 
ton Rouge, Louisiana, USA) Q, AURIGA (Legnaro, 
Italy) g, EXPLORER (CERN) §, NAUTILUS (Fras- 
cati, Italy) (lO) and NIOBE (Perth, Australia) 0. The 
groups operating these detectors agreed in 1997 to start a 
global search for short (~ 1 ms) gw bursts under common 
protocols, by establishing the International Gravitational 
Event Collaboration (IGEG) @. 

All these detectors use the same principles of opera- 
tion. The gw excites the first longitudinal mode of the 
cylindrical bar, which is cooled to cryogenic temperatures 
to reduce the thermal noise and is isolated from seis- 
mic and acoustic disturbances. To measure the strain of 
the bar, a secondary mechanical resonator tuned to the 
cited mode is mounted on one bar face and a sensor mea- 
sures the displacement between the secondary resonator 
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and the bar face. The resulting noise of the detectors in 
terms of strain at the input is 5 — 10 x 10~ 22 /vTIz in a 
bandwidth of ~ 1 Hz surrounding the two coupled-mode 
frequencies. Some of the important physical parameters 
of the five detectors are shown in Table 1 . The detector 
response is optimal for a gw incoming perpendicular to 
the bar axis and polarized along it |0|. The axes of all 
the bar detectors are aligned to within a few degrees of 
one another, so that the chance of coincidence detection 
is maximized. This makes the amplitude acceptance of 
the detectors for the Galactic Center direction greater 
than 0.7 for about 60% of the time Q. 

Each detector output is processed by filters optimized 
for short gw bursts, giving the estimate for the Fourier 
component H(oj) of the strain amplitude h(t) in the de- 
tection bandwidth of ~ 1 Hz around the mode frequen- 
cies listed in Table 1. More specifically, h(t) is the gw 
amplitude multiplied by the antenna pattern of the de- 
tector 1 13 . With the exception of the ALLEGRO detec- 



tor, the noise of the detectors was typically not stationary 
over long observation times and was affected by some un- 
modeled noise sources, whose correlation with common 
environmental noise sources was found to be weak [fL5[ . 
Fig. [j] shows the variability of the noise of each detec- 
tor during 1997-1998, in terms of the Fourier component 
of the gw corresponding to unity signal-to-noise ratio, 
H rms . The detectors had quite similar noise levels, since 
the typical values of H rms were all within a factor of 3. 

We point out that this search for bursts is suitable for 
any transient gw which shows a nearly flat Fourier trans- 
form H(lu) of its amplitude h(t) at the two resonant fre- 
quencies of each detector. The metric perturbation h(t) 
can either be a millisecond pulse, a signal made by a 
few millisecond cycles or a signal sweeping in frequency 
through the detector resonances. The IGEC search is 
therefore sensitive to different kinds of gw sources such 
as a stellar gravitational collapse 0, the last stable or- 
bits of an spiralling N S or BH binary, its merging and 
its final ringdown |16| . The computation of h from the 
measured Fourier component H requires a model for the 
signal shape. A conventionally chosen shape is a pulse 
lasting ~ 10~ 3 s fl7|| ; in this case, H should be multiplied 
by r~j 10 3 Hz to get the corresponding strain amplitude, 
h. 

This letter reports the results of the first coincidence 
search for gw bursts performed by the IGEC observatory. 
The observations covered most of 1997— 1998, including 
625.0 days with at least one detector in operation, 260.4 
days with at least two detectors in simultaneous opera- 
tion, 89.7 days with three detectors, and 15.5 days with 
four. This is the first search with significant observation 
time with more than two detectors. The duration of si- 
multaneous operation would have been greater if it had 
been possible to operate these instruments with higher 
duty factors, which were typically < 50% during this pe- 
riod with the exception of ALLEGRO. More details on 



the observatory, its data exchange protocol and the ex- 
changed data set can be found in Ref. Q . 

The analysis of the data can be divided into two parts: 
a generation of candidate event lists for each of the in- 
dividual detectors, and a time coincidence analysis using 
the lists. This approach, though not optimal, has the 
advantage of being easily implemented and provides for 
a satisfactory effectiveness. 

Each IGEC group extracted the candidates for gw 
bursts, or events, by applying a threshold to the filtered 
output of the detector. The events were described by 
their Fourier magnitude H , their arrival time, the detec- 
tor noise at that time and other auxiliary information. To 
limit the expected rate of accidental coincidences, each 
detector threshold was adaptively set to obtain a maxi- 
mum event rate of ~ 100 / day, with typical values in the 
range Hdet ~2-6x 10~ 21 H z^ 1 corresponding to mag- 
nitude signal-to- noise ratio SNR ~ 3 — 5. Single spurious 
excitations are vetoed against disturbances detected by 
environmental sensors. The AURIGA detector checked 
each event against the expected waveform template by 
means of a x 2 test [(L8) . The lists of the events exchanged 
within IGEC by each detector also include declarations 
of the off- and on- times for the detectors. 

All searches for coincident events used a time window 
of 1.0 second. This choice limits the false dismissal prob- 
ability to less than a few per cent while it ensures a very 
low false alarm probability when at least three detec- 
tors are observing simultaneously. No three- and four- 
fold coincidence was detected, and therefore we did not 
identify candidates for gravitational wave detection in 
the 89.7 days of three-fold observation. The detector 
thresholds were typically 3 x 10 -21 H z~ x for the most 
sensitive three-fold configuration (ALLEGRO-AURIGA- 
NAUTILUS) and 5 x 10~ 21 Hz- 1 for the others. To give 
examples of detectable signals, these thresholds would 
correspond to respectively ~ 0.04 and 0.11 Mq converted 
to isotropic radiation in the optimal polarization at the 
distance of the Galactic Center (10 kPc), assuming a g.w. 
burst of 1 ms duration Jl7| . For comparison, the signal 
expected from the last stable orbits of an optimally ori- 
ented NS coalescing binary at 10 kPc with 2 x 1.4 Mq, 
would give H(lo) = 3 - 4 x lO" 21 ^" 1 at the detector 
resonant frequencies. The number and amplitude of the 
two-fold coincidences found in the 260.4 days of two-fold 
observation are in agreement with the estimated acciden- 
tal background jli] . 

The estimation of the false alarm rate is a crucial el- 
ement in any gw search. It allows for the interpretation 
of any observed coincidences as well as the evaluation of 
the potential of the observatory. Since the events arrival 
times of each detector are randomly distributed with a 
non stationary rate, the expected background of acciden- 
tal coincidences can be computed by two methods: i) by 
modeling the event times as Poisson point process and us- 
ing the measured rates of events for each detector, and ii) 
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by counting the coincidences after performing even time 
shifts of the event times of one detector with respect to 
the others |{|. 

In the first approach, the expected rate of accidental 
coincidences is Ectl 
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where N is the number of detectors simultaneously op- 
erating, T bs their common observation time, At = 1 s 
the maximum time separation for a coincidence, the 
number of events of the ith detector during T b s . This 
equation holds even if the event rates of detectors are 
not stationary as long as they are uncorrelated among 
different detectors. 

The second approach is more empirical. In the case of 
the two-fold coincidence searches, the time shift res ults 
are in agreement with those predicted through Eq. 0.1 
tujl, and demonstrate that the event rates of different 
detectors are uncorrelated. 

The capabilities of the IGEC observatory with respect 
to the false alarm probability are shown in Fig. |^ for a few 
sample configurations of the observatory. The acciden- 
tal rate is calculated as a function of a signal amplitude 
threshold Hthr at the detectors by applying Eq. |0.l| to 
the number of events of the detectors whose amplitude 
is > H t hr- The typical time variability of the instanta- 
neous accidental rate A has been calculated by means of 
a Monte Carlo simulation based on the measured non- 
stationary behavior of event rates on single detectors. 
This variability is about one order of magnitude with 
respect to the mean and is mainly determined by the 
non stationary performances of the detectors. The esti- 
mated mean background of two-fold coincidences is still 
fairly high, unless Hthr is raised well above the data ex- 
change threshold Hdet- On the other hand, a three- fold 
or four-fold coincidence search keeps a high statistical 
significance even for H t hr ~ Hdet-, since the expected 
accidental rates are low enough: respectively less than 
1 false alarm per 10 4 or 10 6 years of observation at 
Hthr ~ 4 x 10 _21 i/z _1 . These fall rapidly as H thr in- 
creases. In fact, the IGEC accidental background noise 
would remain negligible even after centuries of observa- 
tion time. 

The 260 days of observation with two or more detectors 
in simultaneous operation improved by about a factor of 
three the previously set upper limit on the rate of gw 
bursts incident on the Earth S. Assuming the emission 
is described by a stationary Poisson point process and 
using the same procedure as in Ref. , the limiting rate 
set for a gw burst emission from isotropically distributed 
sources is < 4 year^ 1 for H gw > 10~ 20 Hz~ 1 (Fig. ||) with 
95% confidence. A more complete analysis is in progress. 

The IGEC observatory can also be used to set an up- 
per limit on the amplitude of gw bursts corresponding to 



astronomical events, such as supernovae or gamma ray 
bursts. For time windows of the order of the hour or 
larger, each detector is likely to show accidental events 
and therefore this upper limit benefits from a multiple 
coincidence search among the operating detectors. 

A sample of the upper limits on the amplitude of gw 
bursts occurring within a time span of 1 hour, is shown in 
Fig. H for a few weeks of 1998, when up to four detectors 
were operating. These limits apply to the component of 
the radiation emitted with optimal polarization from a 
source optimally oriented with respect to the detectors. 
We are 95% confident that there was no radiation above 
this level hour by hour. In all of 1998, the limits set by 



the IGEC observatory were better than H„ 



6 and 



4 x 10 _21 i/z _1 for 94% and 21% of the year respectively. 
To specialize these upper limits for a specific source direc- 
tion, each detector response should be divided by its an- 
tenna pattern |p| . The corresponding observation times 
of the Galactic Center by IGEC within the same lim- 
its have been respectively 44% and 7.5% of 1998. For 
a source at the Galactic Center emitting isotropically a 
1 ms burst Jj"7| , the above upper limits correspond to 
about 0.16 and 0.07 M & converted in the optimal polar- 
ization. 

Finally, we remark that the IGEC observatory is ca- 
pable of monitoring the strongest galactic sources with 
a very low false alarm probability when at least three 
detectors are simultaneously operating. All the groups 
involved are actively working for upgrading the current 
detector performances and therefore we expect in the 
near future to extend the observation range to the Local 
Group of galaxies, which means an increase of a factor of 
10 of the observed mass. 
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FIG. 1. Spread of the noise of detectors during 1997-1998 
in terms of the Fourier component H rms of the gw at 
SNR = 1. The plotted bands of variability of the noise are 
delimited by selected percentiles, i.e. by selected fractions of 
the observation time for which the sensitivity has been bet- 
ter than H rm3 : bold tick 50%, gray band 16 — 84%, white 
band 2.5 - 97.5%, "T" lines - 100%. The corresponding gw 
amplitude h rms for a ~ 10~ 3 s burst is sketched in the upper 
scale. 



TABLE I. Main characteristics of the IGEC cryogenic bar detectors. The detectors measure the me an Fourier component 
H of the gw in the detection bandwidth of ~ 1 Hz around the mode frequencies. H = (ALu 2 ) -1 *J E / M , where E is the energy 
deposited in the bar by the gw and v is the mean of the mode frequencies. The bars are made by A15056 except for NIOBE, 
whose bar is made of Nb. The sub-kelvin detectors and NIOBE showed very similar typical energy sensitivity in 1997-1998, 
better of a factor of about 4 with respect to the other detectors. The differences in mass and material, though, affect the gw 
sensitivity and give a conversion factor from \fE to H which is 2.3 times worse for NIOBE than for the other detectors. 



detector 


ALLEGRO 


AURIGA 


EXPLORER 


NAUTILUS 


NIOBE 


Mode frequencies [Hz] 


895, 920 


912, 930 


905, 921 


908, 924 


694, 713 


Bar mass M [kg] 


2296 


2230 


2270 


2260 


1500 


Bar length L [m] 


3.0 


2.9 


3.0 


3.0 


2.75 


Bar temperature [K] 


4.2 


0.2 


2.6 


0.1 


5.0 


Longitude 


91°10'44"W 


11°56'54"£ 


Q°12'E 


12°40'21"£; 


115°49'£ 


Latitude 


30° 27' 45" N 


45° 21' 12" N 


46° 27' A 


41°49'26"A 


31°56'5 


Azimuth 


40° W 


44° £ 


39° E 


44° E 


0° 



4 



(D 

-I— ' 



s 1 

io- 2 
10- 4 

10" 6 
10" 8 



10 



10 



10 



12 



m * 














///////////////////A 
















— \- s 


1 
























ALl 


.EG 


RC 


)-AURIGA 










89 d 


ays 


ov 


er a total of 1 04 days; 
















! 
















1 
















! 








s A 


LLE 


3R 


D-f 


OJRIGA-NAUTILUS 5 
















14.3 of 17.1 days I 


















5 ALLE( 


3RCM 


URIC 


■A 










-NAUTI 


LUS-E 


XPLC 


)REF 


\ 








= 6.7 


of 7.9 ( 


Jays 





























3-10 



IO" 20 H thr (Hz" 1 ) 

Estimated rate of accidental coincidences, 



FIG. 2. 

A [year -1 ], versus the threshold H t hr [Hz' 1 ] for a sample 
pair, triple and four-tuple of detectors in 1997- 1998. The 
continuous lines show the mean value of A for signal ampli- 
tudes > Hthr- The dashed lines represent the one std. dev. 
upper bounds for the time variation of the instantaneous ac- 
cidental rates. This figure takes into account the best 85% of 
common observation times, when every detector had an event 
search threshold lower than 3.25, 3.8 and 6.5 x 10 -21 Hz -1 , 
respectively for the pair, the triple and the four-tuple. The A 
for the other operative configurations of detectors were sim- 
ilar, allowing for a small increase of the corresponding H t hr, 
at most by a factor of 2. The bold horizontal line with arrow 
stands for the new upper limit set by all IGEC detectors on 
the rate of incoming gw bursts during 1997-1998. 
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FIG. 3. A sample of the upper limit with 95% confidence on the amplitude of single gw bursts incident with optimal 
polarization and orientation on the IGEC observatory hour by hour in June 1998. The highest peaks shown are due to single 
high amplitude events of one detector while the others were not operating. 
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